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Pohl et al. have reported a very precise measurement of the Lamb-shift in muonic Hydrogen pQ, 
from which they infer the radius characterizing the proton's charge distribution. The result is 5 
standard deviations away from the one of the CODATA compilation of physical constants. This 
has been interpreted [Ij as possibly requiring a 4.9 standard-deviation modification of the Rydberg 
constant, to a new value that would be precise to 3.3 parts in 10 13 , as well as putative evidence for 
physics beyond the standard model [2]. I demonstrate that these options are unsubstantiated. 

PACS numbers: 31.30.jr, 12.20.-m, 32.30.-r, 21. 10. Ft 



I. INTRODUCTION 

The issue is extremely simple. The discrepancy quoted 
in the abstract is between results which do not depend on 
a specific model of the proton's form factor and results, 
by Pohl et al., which do [3]. The conclusion is not that 
the experiments or the theory are wrong, but that the 
model (the customary dipole form factor) is inadequate 
at the level of precision demanded by the data. The 
experiments and QED are right, the dipole is wrong. 

The conclusion of the previous paragraph is the ex- 
pected one. The dipole form factor is but a rough de- 
scription of higher-energy data and is unacceptable on 
grounds of the analyticity requirements stemming from 
causality and the locality of fundamental interactions. 

Moreover, any simple one-parameter description of the 
proton's non-relativistic Sacks form factor, Ge{— q 2 ) in 
terms of only one mass parameter is inaccurate: the pro- 
ton is not so simple. More precisely, the proton's rela- 
tivistic form factor, Ge(q 2 ), is expected, in the timelike 
domain q 2 > 0, to have a complex structure, with a first 
cut starting at q 2 = 4 m 2 and a plethora of branch cuts 
and complex resonant poles thereafter jl]. 

The same is true of the charge distribution, p p (r), the 
Fourier transform of Ge(— q 2 )- Even most naively, p p {r) 
is expected to have a "core" and a "pion cloud" [5], cor- 
responding to a minimum of two length parameters. 



The charge distribution is related to the non- 
relativistic limit of the electric form- factor, Ge, by the 
Fourier transformation 
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Precise measurements of (r 2 ) have two origins. One 
is mainly based on the theory [5] and observations [7] of 
Hydrogen. The result, compiled in CODATA jS], is 



r 2 ) (CODATA) = (0.8768 ± 0.0069 f) 2 



(3) 



The second type of measurement is based on the the- 
ory and observations [3 HH] of very low-energy electron- 
proton scattering. It yields 



(r 2 ){ep) = (0.895 ±0.018 f) 2 



(4) 



This result requires a sophisticated data analysis, partly 
based on a continued- fraction expansion of Ge [5]. 

The two quoted methods of measuring (r 2 ) are model- 
independent, in the sense of not assuming a particular 
form of the proton's charge distribution, p p (r). 

The plot thickens as one considers the Lamb shift 



opF=2 
Zr 3/2 



2Sf /o 1 in the pip atom, measured [Tj to be 
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206.2949 ± 0.0032 meV. 



(5) 



II. IN DETAIL 



In meV units for energy and fermi units for the radii, the 
predicted value [TT] is of the form 



Let £ denote an electron or a p . The leading proton- 
size correction to the energy levels of an ip atom is 
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where (r 2 ) is the mean square radius of p p (r). 
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209.9779(49) - 5.2262 (r 2 



0.00913 (r|) (2) (6) 



The first two coefficients are best estimates of many con- 
tributions while the third stems from the n = 2 value of 
an addend [T2] 
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proportional to the third Zemach moment 
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For a single-parameter description of the charge distri- 
bution, there is an explicit relation between (fp)(2) and 
(r 2 ). Consider, as an example, a p-dominated form factor 
in its narrow-width non-rclativistic limit 
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The corresponding charge distribution is a Yukawian 
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p(r) = 



Its relevant moments are 
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(r 2 )=6/m 2 , (r 3 )=24/m 3 , ( r 3 ) (2) =60/m 3 (11) 



The model-dependent relation is thus 
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For a dipole form factor 
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(q 2 + m 2 ,) 2 
the charge distribution is an exponential 
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for which (r° 

(r 2 ) = 12/m 2 , (r 3 )=60/m 3 , (r 3 ) (2) =315/(2 mg) (15) 
The model-dependent relation is thus 
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The ratio of the numbers in Eqs. ( |12|16[ ) is 441/512 ~ 
0.86, showing the difference of relevant moments between 
to two form-factor "models" . Even if we took the sixth 
root of this number to bring it closer to unity -as exper- 
imentalists do with (r 2 ) to halve the relative error- the 
result would, at the required great precision, still epito- 
mize the model-dependence of the results. 



A. A toy model 



Ge and Gm of protons and neutrons in terms of dis- 
persion relations for the photon propagator involves, lit- 
erally, dozens of parameters [3]- The form- factor "toy 
model" I am going to discuss is not intended to compete 
in accuracy with the dispersive approaches, but only to 
elucidate the current discussion. 

In [1], an accurate description of the theoretically- 
calculated 2ir continuum required products of up to three 
poles. I parametrize p(r) as an interpolation between the 
charge densities of a "p" single pole and a "2-7r" dipole: 
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D = M 2 cos 2 (6)+m 2 sm 2 (0) 
whose two first relevant moments are (r ) = 1 and 
2 6 12 
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To introduce the third Zemach moment, let s = sin(#) 
and c = cos(0). Then 
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We can now check the compatibility of the CODATA 
result of Eq. ^ with the Lamb shift result of Eq. ^ in 
the following way. Solve the two equations 



(r 2 ) (CODATA) = (r 2 
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in M and m for fixed mixing (fixed s 2 ). The results are 
shown in the top Fig. [l] while those of a similar exercise 



with (r 2 ) (CODATA) substituted by (r 2 )(ep) in Eq. ([20] 
are shown in the bottom panel 

From these figures we can draw three conclusions 



1) The system of Eqs. ( |20|2l| ) is soluble only for s 2 >0.1. 
A single-pole or single-dipole are excluded, as expected. 

2) The extracted M and rn are not unreasonable. M 
turns out to be of 0(m p ), while m, which corresponds to 
a dipole paramctrization of "everything but the p pole" is 
not a good enough simplification, a result with m > 2 
would have been nicer. Yukawa intuited pions in a very 
similar manner, but only one at a time. 

3) All experimental results are compatible. 



The photon propagator in the time-like domain (q 2 > 0) 
has led to considerable revolutions (e.g. the discovery 
and interpretation of the J/^), as well as interesting 
challenges, in particular close to its cut at q 2 > Arn^. 
The modeling of the electric and magnetic form factors 



III. CONCLUSIONS 

We face a choice between the following conclusions: 
• The experimental results are not right. 
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FIG. 1: Parameters M and m for which the toy model is 
compatible with the data, with s 2 — sin 2 (8) varying along the 
curves, see Eqs. ( 20|21 1. Top: Lyman in the [ip atom and 
CODATA, shown for the central value and a very asymmet- 
ric ±3(J. Bottom: CODATA substituted for ep scattering, 
central value and ±1 a (there is no solution for + 3 a). 



• The relevant QED calculations are incorrect. 

• There is, at extremely low energies and at the level 
of accuracy of the £p-ad,om experiments, "physics 
beyond the standard model" . 



The combination of the very impressive results in 
Eqs. ( 3|5|6 ) yields a value: 
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3.32 ± 0.22 f 
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r"). The value in Eq. (22 1 looks incredibly large at 



but it is not so unexpected 



with the error dominated by the CODATA uncertainty 
on (- 2 
first 

The result Eq. (22) is p p (r)-independent; to be treated 
with due respect. Right after offering excuses, I shall 
break this rule. The third Zemach moment is very sensi- 
tive to the long-distance part of p(r), compare it to the r 3 
moments of Eqs. ( 11|15 ). Suppose that p{r) has a "core" 
and a "tail" contributing 50-50 to the proton's charge, 
and that the tail's Ge{Ol) is described by a dipole. To 
what scale, m, to does this tail correspond? The value of 
(r 3 }( 2 ) is 1/2 of the one in Eq. (|l5|). Equate it to Eq. (|22j) 



to obtain m ~ 245 MeV, tantahzmgly close to the thresh- 
old of the proton form factor's cut at 2 m n ± ~ 278 MeV. 



r ~ 3.2 f 




q (fm- 1 ) 



FIG. 2: Low-|q| data, compiled and analized by Sick The 
lines are my addition. Only r~0.9 f is bracketed by the data, 
which are all to one side of r ~ 3.2 f. The continuous curve 
-drawn assuming that the absolute data normalization is not 
sacred- illustrates a possible shape whose corresponding p v (r) 



would have a conventional (r p ) and a "large" 
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• A single-dipole form factor is not adequate to the 
analysis of precise low-energy data. 

I have argued that the last choice is the most compelling. 

The theoretical and experimental results I have quoted 
momentarily culminate 125 years of progress in the un- 
derstanding of Hydrogen and its muonic sibling (I am 
setting t = at the discovery date [13] of his famous 
"series" by the Swiss physicist Johann Jakob Balmer). 



IV. DISCUSSION 

Very soon after "v2" of this paper appeared in arXiv, 
a preprint by Cloet and Miller was posted [T3] . These 
authors find it impossible to extract a result as large as 
that of Eq. (22 1 from ep scattering data. 



A crucial problem in this connection was adroitly em- 
phasized by Sick 9 . It is very difficult to extract reliable 
information on p(r) from its Fourier transform, G^;(q) 2 . 



4 



The radius of convergence of the expansion in r from 
which one extracts (r^) is so small, that one must use sim- 
ulations and a continued-fraction expansion to obtain a 
stable, numerically-meaningful result not contaminated, 
for instance, by the term in (rz). Clearly, if extricating 



(r 2 ) is delicate, the more so it is to infer 



(^>(2)- 



A problem in extracting (^p)(2) from ep data is illus- 
trated in Fig. [2j borrowed from [5]. The data amply 
bracket a domain around r ~ 0.9 f, required to measure 
an 



7-p) 1 / 2 of this order. Contrariwise, these data are all 



to one side of 3.2 f, the (r;)^ scale in Eq. (|22|. A 

data 



result based on them has to be an extrapolation 



of 



with a large spread and a poor \ P er degree of freedom. 

The absolute normalization of the data at small q 2 
(and their always disdained systematic errors) are a no- 
torious hurdle [9j. If one tolerates a few % uncertainty, 
the actual Gg(q) may correspond to the shape shown in 
Fig. [2| reminiscent of the tttt contribution to F1.2 calcu- 
latedin [J. The Fourier transform of this shape would 



have the customary (r 2 ) but a "surprisingly" large (r|)(2). 

Theoretical estimates of (»"p)(2) would be interesting. 
In their study [I], Belushkin et al. find that the 2-k con- 
tinuum is a very significant contribution, rising dramati- 
cally below its scale, |q| = 2 ~ 0.71 f -1 , and requiring 
for its parametrization products of up to three poles. 
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